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Baryon stopping at the SPS and RHIC energies is calculated by introducing a new 
baryon junction mechanism into HIJING. The exchange of a baryon junction, according 
to Regge phenomenology, leads to a cosh[(y — ?/cm)/2] rapidity dependence and an 
energy dependence of the inclusive baryon cross section. This baryon junction dynamics 
also leads naturally to enhanced pr broadening in pA and AA together with enhanced 
mid-rapidity hyperon production. 
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1. Introduction 

The systematic study of baryon stopping and hyperon production in pp,pA, and AA 
collisions is essential in order to differentiate the sought after new physics of dense equili- 
brated matter from non-equilibrium multi-particle production dynamics [p]]. Baryon stop- 
ping refers to the transport of baryon number in rapidity space away from the nuclear 
fragmentation regions and is measured through the single inclusive rapidity distribution 
of protons and hyperons. As discussed in fHH, recent SPS data shows a high degree of 
baryon stopping and anomalous hyperon production. These findings are based on com- 
parisons of predictions using the HIJING || and VENUS Monte Carlo event generators 
with SPS data for p + A and A + B collisions. In particular, the baryon stopping mech- 
anism modeled in HIJING was found to underestimate the stopping in heavy nuclear 
collisions, while VENUS model was able to reproduce the baryon stopping and strange 
baryon production by introducing a new non-equilibrium dynamical mechanism called 
the double string. The mechanism in HIJING, adapted from the LUND Fritiof J7| and 
dual parton model (DPM) ||, involves diquark-quark string excitations followed by string 
breaking. The final valence baryons in this model always emerge within about a unit of 
rapidity from those of the unbroken diquarks. This association of the baryon with the 
end point diquark creates a characteristic baryon rapidity distribution where the baryons 
are concentrated near the beam and target fragmentation regions as shown in Fig la. 
The extrapolation of the pp physics of HIJING to PbPb is seen to underestimate the 
mid-rapidity valence baryon yield by a factor of two in Fig lb. VENUS avoids this "fly 
trap" by postulating a specific model of diquark breakup, leading to double strings. The 
double string emerging from a broken diquark is assumed to loop around a valence quark 
from the other fragmentation region. The produced baryon is then associated with the 
fragmentation of the double string near its 180° bend. The parameters of this mechanism 
are adjusted from pA and suffice to reproduce the AA data. We emphasize that the 
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analytic one parameter Multi-Chain- Model |4| also reproduces well [[3] the valence proton 
distributions based on fitting pA. The advantage of a Monte Carlo event generator is the 
ability to test via other observables such as hyperon production and p± systematics the 
consistency of specific dynamical assumptions with a wide range of data. 

Another model of baryon transport was proposed by Kopeliovich and Zakharov || based 
upon pQCD considerations. That mechanism has been recently included in a new version 



of the DPM [I0|. In that model, the leading diquark is broken by a color exchange which 
changes its color state from a {3} to a {6}. This allows one of the valence quarks to flow 
into the central rapidity region and act as the seed for the production of the final baryon 
number. The inclusive cross section of this mechanism is characterized by a coshy/2 
rapidity dependence and an energy dependence which follows from the assumed 

1/ \fx distribution of the valence quarks in a proton. In practice this mechanism is similar 
to that invoked in VENUS in that the fragmentation of the now two independent quarks 
from the diquark is similar to that from two independent strings. As in VENUS this 
configuration allows for a factor of 2 enhancement in the production of strange baryons. 
Adding a component of about 20% for each nucleon-nucleon collision, with the assumption 
that once the diquark is broken then it remains broken, DPM is then also able to reproduce 
the SPS data. 



In this work we test a third model of baryon transport recently proposed by Kharzeev|TT 
and based on a long dormant Rossi- Veneziano baryon junction Regge exchange model [|T^ . 
In this approach, the baryon number is traced by a non-perturbative baryon junction of 
topological nature. The special theoretical appeal of that idea arises from the form of 
the baryon wavefunction in QCD. The requirement of gauge invariance of the nonlocal 
operator creating a baryon naturally leads to the concept of a baryon junction: 



B(xi,x 2 ,x 3 ,xj) = e ijk 



Pexp^ig J dx^A^jq(xi) Pexp^ig J dx^A^j g(x 2 

P exp \ig J dx^A^j q(xs) 



• (1) 

k 



The baryon junction is a vertex at xj where the three gluon Wilson lines emanating from 
the three valence quarks (in SU(3)) must join in order to form a gauge invariant operator. 
In a highly excited baryonic state, the three valence quarks fragment via multiple qq into 
mesons leaving three sea quarks eventually around the junction to form the observed 
final baryon. This is the sense in which the junction traces the baryon number. Being a 
purely gluonic configuration, the junction may be easily transported into the mid-rapidity 
region. The fragmentation of the end point valence quarks naturally leads to three beam 
jets, very similar to the previous mechanisms considered. The primary advantage of this 
mechanism is that Regge phenomenology can be used to estimate the inclusive inelastic 
cross section for this process. From pp and pp data, the junction-anti-junction exchange 
is characterized by a trajectory with intercept, ajj(0) ~ 1/2. This leads immediately to 
a coshy/2 rapidity dependence and an energy dependence. 

The unique prediction of this mechanism is that fragmentation of the valence quarks 
down to the junction is expected to enhance hyperon production by a factor of 3 just 
from the random combinatorics of s vs u, d pair production. In addition, from the random 
addition of three sea quarks, the transverse momentum of the final baryon is automatically 
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Figure 1. A comparison of HIJING (solid), HIJING/B (dashed) and HIJING/B with 
"ropes" (dotted) with various experimental data. 



enhanced by a factor of y/3. 

Given these interesting coupled consequences of the baryon exchange mechanism, we 
have chosen to implement this model in a new version of HIJING, called HIJING/B ITBl. 



2. HIJING/B 

In HIJING/B, baryon junction stopping is implemented using a "Y" string configuration 
for the excited baryons in which the baryon junction and produced baryon are distributed 
as noted above. The resulting three beam jets are treated as qq configurations. For each 
nucleon-nucleon interaction, a ~ 60% probability is given that only one of the two nucleons 
will have its baryon junction stopped in the central rapidity region. Baryon excitations 
without junction exchange are taken to be standard (qq) — q strings for fragmentation 
purposes. Due to kinematical considerations, junction exchange is allowed only if the 
invariant mass exceeds m > 4.5 GeV so that the three beam jets have room to decay. At 
SPS energies this trivial kinematic constraint limits considerably the number of junction 
exchanges allowed. The junction exchange probability is fitted so as to reproduce the pp 
valence proton data. For multiple collisions, we assume that once the baryon junction is 
exchanged, it remains stopped in subsequent soft interactions. 

A comparison of the prediction of HIJING and HIJING/B with SPS data [|l4l-4l~6t1 is 
given in Fig la,b,c. In part (a), (b) we show the valence proton (p—p) rapidity distribution 
for, S+S and Pb+Pb, respectively. Whereas HIJING under-predicts the baryon stopping, 
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HIJING/B is shown to provide a large amount of stopping. The valence hyperon (A — A) 
rapidity distribution is shown in part (c) of Fig. 1. Here, HIJING/B provides a factor of 
3 more valence A at mid-rapidity consistent with the expected strangeness enhancement 
associated with the baryon junction exchange. However the very large value of the valence 
hyperon yield in PbPb cannot be reproduced. Multiple final state interactions or "rope" 
effects must be included to understand those yields. Increasing the effective string tension 
by a factor of two makes little difference on the valence non-strange baryon yields, but 
can account for the anomalous hyperon yield (dotted curves). 

Finally in Fig Id, the predictions of this model for the valence proton rapidity distri- 
bution in Au + Au collisions at RHIC energies (y/s = 200 GeV) are shown. HIJING/B 
predicts nearly twice the initial baryon density, fm) « 2po, at mid-rapidity than 

HIJING. 

In conclusion, a novel baryon junction stopping mechanism is implemented in HI- 
JING/B to reproduce baryon transport in pp collisions. This is sufficient to account for 
the midrapidity valence baryon data in pA and AA. The model naturally leads to large 
enhancement of the p± slope as well as of the yield of hyperons. However, the preliminary 
A — A yield in PbPb is still underestimated without multiple final state interactions. 
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